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Preface 

A  need  for  atlas -type  climatology  has  prompted  a  re-examination  of  the  map 
projections  that  have  been  used  for  hemispheric  and  world-wide  climatic  presen¬ 
tations.  At  the  same  time  there  has  been  a  need  for  mesoscale  regional  studies  of 
climatology,  which  requirement,  for  this  report,  has  been  viewed  as  a  necessity 
for  improved  geographic  presentations.  Both  types  have  been  met  with  the  develop¬ 
ment  of  equal-area  mapping. 

In  this  effort  we  have  been  assisted  by  the  principal  investigator,  Mr.  Jack 
Mettauer,  under  a  Regis  College  contract,  with  computer  solutions  and  plotting 
of  regional  equal-area  maps. 
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Mapping  the  Climate 


I.  INTRODUCTION 

The  main  thrust  of  this  paper  is  to  detail  a  graphical  means  of  presenting  geo¬ 
graphic  information  on  climate  for  expeditious  use.  Climate,  as  defined  by 
Landsberg^  is  the  "collective  state  of  the  earth's  atmosphere  for  a  given  place 
within  a  specified  interval  of  time."  The  "place,"  as  Landsberg  explains,  may 
be  as  small  as  an  ant  hill,  as  large  as  a  continent,  or  in  fact  the  whole  earth. 

This  paper  is  conceived  as  a  step  toward  the  development  of  mesoscale 
climatology,  or  the  study  of  the  climate's  temporal  and  spatial  variability,  within 
areas  from  a  few  square  kilometers  to  several  thousand  square  kilometers. 

Studies  of  mesoscale  climatology  are  numerous  and  varied  partly  because  they 
need  much  adaptation  to  time  and  place.  In  the  opinion  of  this  writer,  they  do  not 
permit  practical  generalizations  for  routine  use.  As  a  consequence,  each  region 
small  or  large,  requires  an  individual  description  of  its  climatic  characteristics. 
Yet  for  operational  and  planning  purposes,  it  is  important  that  the  climate  be 
presented  in  sufficient  detail  to  distinguish  important  spatial  differences,  such  as 
the  increase  in  the  frequency  of  cloud  cover  from  the  coastline  toward  the  foothills 
of  a  west  coastal  region.  An  example  of  the  effects  of  a  large  body  of  water,  and 

(Deceived  for  publication  l(i  January  1981) 

1.  Landsberg,  Helmut  (1958)  Physical  Climatology,  Gray  Printing  Co.  ,  Inc., 
DuBois,  Pennsylvania,  44fl  pp. 
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the  distance  of  the  station  therefrom,  can  be  perceived  on  diurnal  wind  variations 
2 

east  of  Lake  Michigan. 

The  practice  of  presenting  climatic  information  in  map  form  is  in  line  with 
r-’i’ent  developments  in  computer  graphics.  A  compilation  of  papers,  consisting 

"of  ideas,  practical  methods  and  procedures,  and  illustrated  applications,"  has 

3 

been  published  recently  by  the  Harvard  Libra rv  of  Computer  Graphics.  Among 

4 

these  papers,  one  by  Box  describes  several  ecological  studies,  which  are  facil¬ 
itated  by  a  system  of  " sytlagraphic"  mapping,  anachmnized  as  SYMAP.  Among 
the  environmental  elements  presented  graphically  are  climatic  factors  including 
temperature,  precipitation,  evaporation  and  solar  radiation.  Box  describes  the 
use  of  a  new  projection  by  Robinson0  which  represents  a  compromise  between  equal 
areas  and  correct  shapes. 

Programmed  procedures  for  the  computer-plotting  of  boundaries  and  climatic 
information  onto  a  map  have  been  described  by  others,  with  systems  generally 
given  code  names.  Recently  a  system,  called  EASYMAP  has  been  developed  to 
depict  statistical  information  by  geographic  areas,  using  shading,  technically  called 
choropleth.  Present  users  include  the  Bureau  of  the  Census  and  the  United  States 
Dept,  of  Commerce.  While  more  attractive  systems  may  be  available  for  such 
mapping,  aesthetics  is  sacrificed  for  economy  and  availability.  The  boundaries  of 
geographic  entities  such  as  states  or  countries  are  plotted  by  the  coordinates  of 
many  points,  before  the  climatic  information  is  plotted.  Data  sets,  for  this  purpose, 
for  the  automatic  plotting  of  shorelines,  rivers,  lakes,  canals,  political  boundaries 

7 

are  now  available. 

Samples  of  mesoscale  mapping  appear  often  in  t he  meteorological  literature. 

Q 

In  Swan  and  Lee,  mapping  is  done  in  a  rectangular  area,  including  San  Francisco, 


2.  Weber,  M.R.  (1978)  Average  Diurnal  Wind  Variation  in  Southwestern  Lower 

Michigan,  J.  Appl.  Meteorol.  17:1182- 1189. 

3.  Harvard  University  (1980)  The  Harvard  Library  of  Computer  Graphics  Mapping 

Collection.  Vol.  1—11,  Center  for  Management  Research,  850  Boylston  St., 
Chestnut  Hill,  Massachusetts  02107. 

4.  Box,  Elgene  (1980)  Use  of  Svnagraphic  Computer  Mapping  in  Geoecology,  The 

Harvard  Library  of  Computer  Graphics.  1979  Mapping  Collection,  Vol.  a, 
pp  11-27. 

5.  Robinson,  A.H.  (1974)  A  New  Map  Projection:  Its  Development  and  Charac¬ 

teristics,  International  Yearbook  of  Cartography  14:145-55. 

(!.  Iskow,  L.  I.  ,  and  Spoeri,  R.K.  (1979)  EASYMAP:  A  Machine-Independent 
Line  Printer  Statistical  Mapping  System,  The  American  Statistician 
33:223-224. 

7.  Herman,  Alan  (1980)  Computer  Mapping  Applications  at  the  Atlantic  Oceano¬ 

graphic  and  Meteorological  Laboratories,  Harvard  Library  of  Computer 
Graphics,  1980  Mapping  Collection,  Vol.  8^  pp  105-181. 

8.  Swan,  P.  R.  ,  and  Lee,  I.Y.  (1980)  Meteorological  and  Air  Pollution  Modeling 

for  an  Urban  Airport,  J.  Appl.  Meteorol.  19:534-544. 


Ill  km  north-south  and  88  km  east -west.  In  this  small  area  are  measurements 

made  at  some  13  sites,  resulting  in  "interesting"  mesoscale  patterns  of  winds  and 

air  pollution.  In  addition,  macroscale  computer  mapping  has  been  illustrated  with 

g 

upper-air  data  in  the  southern  hemisphere. 

In  mapping  we  must  decide  on  the  nature  of  the  statistical  entries,  in  order  to 
provide  graphically  information  on  the  climate.  The  most  obvious  entry  is  a  num¬ 
ber  properly  located  on  the  map.  However,  the  purpose,  the  providing  of  informa¬ 
tion,  may  be  served  better  by  isopleths  or  isolines  of  equal  values  of  the  climatic 
element  or  parameter. 

There  are  depictions  of  statistical  information  by  shading  (choropleth)  or  by 
colors.  An  example  of  shading  can  be  found  in  the  Climatic  Atlas  of  the  United 
States,  ^  published  by  the  United  States  Department  of  Commerce;  for  an  example 
of  the  use  of  colors,  the  Climatic  Atlas  of  Australia.  1  issued  by  the  Australia 
Government  Publishing  Service  may  be  cited.  With  the  development  of  computer 

science  and  equipment,  it  is  now  possible  to  plot  information  automatically  by  use 

12 

of  shading  or  colors,  for  optimum  presentation  and  use. 

There  are  advantages,  but  also  obstacles  in  the  use  of  graphics,  especially 
for  multidimensional  data,  such  as  ceiling  and  visibility;  in  fact,  a  satisfactory 
means  of  measuring  the  utility  of  a  graphical  scheme  is  still  lacking.  A  fairly 

reliable  element  of  personal  judgment  does  exist  in  the  selection  of  shading  or 

13 

gradation  of  color  intensities.  However,  Wainer  and  Franrolini  think  that  a 
legend  of  shading  or  colors  is  a  "venial  sin"  if  it  must  be  memorized  by  the  viewer 
in  order  to  read  the  map. 

2.  ON  THK  CHOICE  OF  MAP  PROJECTION 

14 

All  maps  inevitably  distort  the  shape  of  a  global  land  mass.  The  earth's 
sphere  cannot  be  mapped  accurately  onto  a  flat  surface,  although  the  approximation 

9.  Levanson,  N. ,  and  Julian,  P.  R.  (1980)  Antarctic  laO  mb  Pressure  Maps  from 
TWERI.E  and  Radiosondes  (November  1975-March  197(>),  Monthly  Weather 
Review  108:520-52(>. 

10.  U.S.  Dept,  of  Commerce  ( 1908)  Climatic  Atlas  of  the  United  States.  ESSA 

EDS,  U.S.  Government  Printing  Office,  Washington,  D.  C. 

11.  Bureau  of  Meteorology  (1977)  Climatic  Atlas  of  Australia.  Australia 

Government  Publishing  Service,  Canberra. 

12.  Feinberg,  S.  E.  ( 1979)  Graphical  Methods  in  Statistics,  The  American 

Stat ist ician  33;lt>5-178. 

13.  Wainer.  I).,  and  Franrolini,  R.J.  (1980)  An  Empirical  Inquiry  Concerning 

Human  Understanding  of  Two-Variable  Color  Maps,  The  American 
Statistician  3jl:8 1  -93 . 

14.  Raisz,  E.  (1902)  Principles  of  Cartography.  McGraw-Hill,  New  York. 

315  pp.  ~ 
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can  bo  good  enough  in  a  small  area  to  make  distortion  almost  imperceptible. 
Virtually,  all  global  maps  must  slimt  the  earth's  surface  spl:’,  usually  along  one 
or  more  meridians.  There  are  two  aspects  to  avoid,  distortion  and  splits,  which 
themselves  may  dictate  the  choice  of  a  map.  Ilut  there  are  other  factors  to  be 
dictated  by  the  map's  purpose. 

Maps  should  have  one  or  more  features  designat'd  ns  follows:  (1)  azimuthal, 
(21  equidistant,  (3)  equal-area,  and  (4)  conformal.  I  di.-cus.-  tln-se  briefly: 

1.  In  an  azimuthal  map,  there  is  a  central  point  (F  ),  :iot  necessarily  repre¬ 
senting  the  north  pole,  although  it  often  does.  From  the  central  point,  in  all  direc¬ 
tions,  points  of  the  same  earth  distance  from  P  will  lie  on  a  circle  (entered  on 

P  Included  among  the  azimuthal  maps  are  the  polar  stereographic  map  and  the 
polar  Lambert  azimuthal  equal-area  map. 

2.  An  equidistant  map  is  also  azimuthal.  The  distance  from  the  central  point 
(P  )  of  the  map  to  another  point  (P)  is  an  exact  representation  of  the  earth  distance 
from  to  P.  This  is  a  good  map  to  have  in  a  ccn'ral  airport,  in  order  to  estimate 
the  true  distance  to  another  airport. 

3.  In  an  equal-area  map,  the  sizes  of  all  areas  are  accurately  portrayed  rela¬ 
tive  to  other  areas  of  the  global  surface.  Such  maps  comprise  a  large  family  of 
projections,  from  which  a  selection  can  ho  made  to  meet  other  rcquiremi  nts  such 
as  economy  of  page  space,  freedom  from  distortion  in  selected  locations,  or 
absence  of  splits  in  important  areas. 

4.  A  conformal  map,  while  conserving  proper  shapes  of  small  regions  such 
as  a  city,  produces  distortion  of  large  areas  such  as  the  continents.  But  because 
of  the  freedom  from  distortion  at  single  points,  conformal  projections  arc  the 
maps  of  choice  in  navigation,  supported  by  gnommiic  charts  for  detailing  the  great 
circles  or  routes  of  shortest  distance. 

For  statistical  purposes,  such  as  ours,  the  equal-area  group  of  maps  offers  a 
distinct  advantage.  At  the  same  time,  there  is  an  advantage  when  true  land  shapes 
are  closely  preserved,  or  when  the  distortion  is  lessened.  Robinson'*  lias  presented 
a  persuasive  case  for  the  general  goal  of  creating  new  map  projections  with  specific 
objectives,  even  though  hundreds  of  projections  now  exist.  Most  cartographers 
have  not  engaged  in  such  efforts  because  of  other  claims  upon  their  time.  However, 
with  the  computer-technology  breakthrough,  it  is  now  much  easier  to  generate  new 
or  novel  base  maps,  and  to  use  them  readily  by  the  same  computer-plotting  tech¬ 
niques. 

The  maps  described  in  this  paper  are  consistently  equal-area.  The  earth  is 
assumed  to  be  absolutely  spherical,  with  exactly  t>0  nm  to  one  degree  of  latitude  or 
great  circle.  To  be  used  in  the  small  areas  of  mesoscale  studies,  the  shapes  of 
azimuthal  charts  arc  correct  and  conformal  at  map  center.  From  the  center 
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outward,  the  distortion  remains  less  than  1  percent  even  as  far  out  as  700  nm  from 
the  center  to  be  discussed  later. 

The  earth's  radius  (R)  lias  value: 

K  -  (180  X  60/ir)  nm  =  3,437.75  nm 

-  3,  958.  02  st  mi 

-  li,  3  (if>.  71  km 

^  2. 50818  X  108  inches 

In  an  equal-area  map,  the  linear  scale  will  vary  with  orientation,  but  the  geometric 
mean  scale  (s)  will  be  constant  throughout  the  map.  Map  sizes  are  given  by  the 
reciprocal  of  the  map  scale,  as  follows: 

Scale 
s 

1:72,  900 
1:729,  (iOO 
1:7,  290,  000 

o  r 

s 

1:100,  000 
1:1,  000,  000 
1:10,  000,  000 

The  symbol  for  latitude  is  c,  positive  for  latitude  north,  negative  for  latitude 
south.  The  symbol  for  longitude  is  A,  positive  when  measured  east  of  the  Greenwich 
meridian,  negative  when  measured  west  of  the  Greenwich  meridian.  The  spherical 
coordinates  (<j>,  A)  are  symbolized  (<>(),  A  for  the  central  point  of  the  azimuthal  equal- 
area  chart. 

The  maps  described  in  this  report  are:  (1)  A  regional  map,  centered  on  a 
selected  point  A  ),  produced  to  meet  the  requirements  of  a  mesoscale  study. 

It  is  always  azimuthal  equal-area.  (2)  World  maps  centered  on  the  north  pole  are 
of  two  types:  First,  a  square  map  that  is  polar  azimuthal  equal-area  for  the  northern 
hemisphere,  but  fits  the  remaining  southern  hemisphere  into  the  completed  square, 


Kilometer  to 
the  Centimeter 

1 

10 

100 


Map  Size 
Nautical  Miles 
to  the  inch 

1 

10 

100 


13 


for  economy  of  space.  Second,  a  world  map  centered  on  the  north  pole,  split  along 
four  meridians  so  as  to  produce  four  quadrants,  three  of  which  display  the  con¬ 
tinents  unsplit,  and  with  relatively  little  distortion  except  for  Antarctica.  The  map 
f'ts  neatly  onto  an  (8-1  2  X  11)  page  of  an  atlas  or  other  publication  of  conventional 
shape. 


AZ1MI  THAI.  Kyi  VI.  ARF.A  MAI’  CKNTKRKD  ON  <b  ,A  \ 

Y  ,vo’  tr 

3.1  (Construction 

In  all  maps  of  this  kind,  there  is  a  central  point  corresponding  to  the  terres¬ 
trial  coordinates  (0^,  A  )t  around  which  the  map  is  drawn,  preserving  the  equal- 
area  feature.  There  is  circular  symmetry  around  the  central  point  (Figure  1)  at 
which  there  is  zero  distortion.  Distortion  gradually  increases  outward  from  the 
central  point.  When  the  projection  is  centered  on  the  north  pole,  it  is  railed 
Polar  Lambert  Azimuthal  Equal-Area. 


Figure  1.  To  Illustrate  the  Cartesian  Coordinates 
of  a  Point  (P)  on  an  Azimuthal  Map  Centered  on  a 
Point  (P  )  Whose  Geographical  Coordinates  are 
($o«  All  points  represented  on  the  circle  are 

geographically  equidistant  from  P  (<J  ,  \  ) 

o  o  o 
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On  a  spherical  earth  (Figure  2)  the  angular  measurement,  in  radians,  of  any 
side  of  a  spherical  triangle  (NPqP)  gives  the  earth  distance  along  the  great  circle 
after  multiplication  by  the  earth's  radius  (R).  The  full  length  of  the  equator,  like 
any  great  circle  is  2nR. 


S 


Figure  2.  To  Illustrate  the  Spherical  Triangle 
NPQP  on  a  Sphere,  in  Which  PD  Has  Coordi¬ 
nates  and  P  Has  Coordinates  (</>,  A.) .  The 

distance  between  P  and  P  is  measured  by  the 
great-circle  arc  (v) 


2 

On  an  equal-area  map  of  scale  (s)  the  area  of  a  global  hemisphere  (2itR  )  is 
mapped  into  the  area  of  the  circle  of  radius  r  such  that 

2?rR2s2  =  trr2 
o 

which  makes 


r  =  v/2  •  (Rs) 
o 


(1] 


Example:  If  the  geometric  mean  scale  (s)  is  1:35,  000,  000  as  it  is  in  the  1070 
working  chart  published  by  the  Defense  Mapping  Agency  (Figure  3)  then  r  is 
10.  135  inches;  that  is,  the  northern  hemisphere  is  portrayed  within  ;<  circle  of 
radius  10.  135  inches. 


Figure  3.  The  Pola.  Equal-Area  Map  of  the  World.  Prepared  at  request  of  the 
Air  Force  Geophysics  Laboratory  by  the  Defense  Mapping  Agency  (PEA -3),  June 


Consider  the  mapping  of  the  global  area  surrounding  the  point  ($Qf  Xq) .  The 
small  circle  that  subtends  the  conical  angle  y  at  the  center  of  the  earth  encloses 
the  area  2ttR“(l  -  cos  y).  Hence  the  radius  (r)  of  the  corresponding  circle  on  the 
map  is  such  that 

2irR“(l  -  cos  v)  •  s^  =  trr^ 


which  makes 


r  =  /v^2(l  -  cos  y)  •  (Rs) 


(2) 


On  the  earth  (Figure  2),  for  the  spherical  triangle  (NPqP),  the  following  equa¬ 
tions  hold: 


sin  A  sin  B 
cos  d  "  cos  <£o 


sin  (X  -  X^> 
sin  y 


(3) 


cos  y  =  sin 


</>q  sin  <j>  +  cos  cos  <j>  cos  (X  -  X^) 


(4) 


sin  <£  =  sin  4(J  cos  y  +  cos  $ sin  y  cos  A  (5) 

On  the  map  (Figure  1)  the  rectangular  Cartesian  coordinates  (x,  y)  correspond¬ 
ing  to  (<t>,  X)  are  given  by 


x  -  r  . 

sin  A 

(fi) 

y  =  r  • 

cos  A 

If  A  >  90°,  then  cos  A  <  0,  to  make  y  negative.  For  A  =  90°,  Eq.  (5)  gives  the 
corresponding  <(> -value  as: 

sin  i  -  sin  <£  cos  y 
s  o  s 

so  that  must  have  the  same  sign  as  d  for  yg  <  90°.  Or  both  <b,  <t>  are  in  the 
same  hemisphere.  Together  with  Kq.  (3) 


cos 


■> 

1  -  sin”  $ 

_ Li. 

.2 

-sin  b  ■  sin 

(  » 


(X  -  X  ) 

o 


(7) 
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For  A  >  90°,  <t>  <  <i>  for  positive  or  negative  <t>  .  Or,  conversely,  when  4>  <  <t>  , 
s  s  s 

then  A  is  taken  as  equal  to  ”  minus  the  acute  angle  given  by  Eq.  (3),  making  the 
y-component  negative. 

Example:  Suppose  we  wish  to  draw  the  parallels  and  meridians  surrounding 
Boston,  Mass.  (42°22'N,  71°2'W).  Table  1  shows  the  results  of  computer  compu¬ 
tation  for  the  coordinates  (x,  v)  given  the  spherical  coordinates  {<t>,  *•).  The  computer 
was  programmed  (Appendix  A)  to  solve  for  (x,  y)  by  the  above  equations  in  the  order: 
Eqs.  (4),  (3),  (7),  (2),  (6). 

Table  1.  Map  Coordinates  Corresponding  to  the  Given  Point  of  Intersection 
of  Meridians  and  Parallels  Surrounding  Boston,  Massachusetts,  as  Com¬ 
puted  with  R  =  2.  508  X  108  inches,  s  =  1:614,  600,  for  <t>  :42°22'N;X  :71°2'W 
1  ’  ’  '  '  ’o  o 


Latitude 

Longitude 

X 

(inches) 

y 

(inches) 

42°00'N 

70°30'w 

2.  82 

-2.  60 

7 1°00'W 

0.  18 

-2.  61 

71°30'W 

-2.  47 

-2.  60 

42°30'N 

70°3 O'W 

2.  80 

0.  95 

7 1°00'W 

0.  18 

0.  95 

7 1**3 1  'YV 

-2.  45 

0.  96 

Table  2  shows  the  computed  Cartesian  coordinates  for  selected  stations  in  the 

Boston  area,  which  are  plotted  in  Figure  4.  Figure  4  retains  the  square  grid  by 

2 

which  all  points  were  plotted.  Each  square  corresponds  to  243.  4  km  of  terrain. 
Another  example  of  Cartesian  coordinates  (x,  y),  found  to  correspond  to  geographic 
coordinates  (<t>,  A)  appears  in  Table  3  and  plotted  in  Figure  5.  It  illustrates  the 
effectiveness  of  the  system  in  equatorial  regions,  around  the  international  date  line, 
in  fact  anywhere.  However,  to  avoid  computational  embarrassments  (such  as  zero 
divided  by  zero)  the  central  point  (<t>o,  A.  )  was  taken  north  of  the  equator  31  meters, 
corresponding  to  1  second  of  angle.  The  plots  in  Figures  6,  7,  8,  9  show  meridians 
and  parallels  that  were  plotted  and  drawn  by  computer  automatically.  The  program 
was  devised  and  written  by  the  analyst,  Jack  Mettauer,  under  a  contract  with 
Regis  College. 
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Station 

d> 

X 

X 

— 

Natick 

42° 17' N 

7 1°2  l'W 

-1.  G7 

-0.  59 

Franklin 

42°5'N 

71°24'W 

-1.  94 

-2.01 

Lexington 

42°27'N 

71°14'W 

-1.05 

0.  59 

Needham 

42°17'N 

71°14'W 

-1.  05 

-0.  59 

Brockton 

42°4'N 

71°  l'W 

0.  09 

-2.  14 

S.  Weymouth 

42°9'N 

70°5fi'W 

0.  53 

-1.  54 

Concord 

42°28'N 

7 1°2  l'W 

-1.  60 

0.  72 

Figure  6  was  drawn  for  Germany,  centered  on  (51“N,  10“E)  w  ith  two  locations 
plotted  (B  for  Berlin,  M  for  Munich).  Figure  7  was  drawn  for  Ecuador,  centered 
on  (1°30'S,  78”00'W),  with  the  locations  of  Lasso  and  Queto  plotted  by  computer. 

Figure  8  was  drawn  for  Australia  in  the  Southern  hemisphere  with  central 
point  at  25°S,  135°30'E.  Figure  9  was  drawn  also  for  Australia,  but  the  azimuthal 
map's  center  was  chosen  at  40°N,  100°E.  Clearly  this  map,  while  still  equal-area, 
is  no  longer  conformal,  even  approximately. 

3.2  Map  Distortion 

On  the  sphere  (Figure  2)  the  incremental  earth  distance  caused  by  an  increase 
(6r  radians)  in  the  arc  (r)  is  11  •  dir.  The  incremental  distance  caused  by  an 
increase  (6 A )  in  the  angle  A  (radians)  is  It  •  sin  r  •  6A.  When  the  two  are  equal, 

6A  -  6r/sin  v  (8) 

On  the  map  (Figure  1)  the  two  corresponding  increments  will  generally  not  be  equal, 
and  their  ratio  will  be  a  measure  of  the  distortion  (e)  thus: 


1  +  e  =  r  •  6A/6r 

which,  together  with  Eq.  (2),  its  differentiation  and  Eq.  (8),  gives 
e  =  2/(1  +  cos  r)  -  1 


Hence,  at  P  where  y  =  0,  there  is  no  distortion.  But  the  distortion  increases  as 
v  increases.  It  will  reach  1  percent  when  y  =  11.42°  which  corresponds  to  an  earth 
distance  of  085.  27  nm,  or  789.  1  st.  miles  or  1,  269.  12  km. 


Table  3.  Cartesian  Coordinates  (x,  y)  for  the  Azimuthal  Equal-Area  Map  of  the 
South  Pacific,  from  10°N  to  10°S,  170°E  and  170°W,  Centered  on  (00°00'01"N, 
180°00'00"  W) 


X 

y 

Name 

Latitude 

Longitude 

inches 

inches 

Gilbert  Island 

00°00' 

174°5'E 

-1.  29 

0.00 

Funafuti  Atoll 

9°00'S 

179°00'E 

-0.  22 

-1.  97 

Nukunono 

9°12'S 

171°55'W 

1.  75 

-2.02 

Canton  Island 

2°46'S 

171°43'W 

1.  81 

-0.  61 

Gardner 

4°41'S 

174°34'W 

1.  19 

-1.03 

Hull 

4°30'S 

172°14'W 

1.  69 

-0.  99 

Beru 

1°21'S 

175°58'E 

-0.  88 

-0.30 

Tarawa 

1°21'N 

172°56'E 

-1.  55 

0.30 

N'anumea 

5°39'S 

176°0R'E 

-0.  85 

-1.  24 

Nui 

7°1G'S 

177°10’E 

-0.  62 

-1.  59 

Baker  Island 

00°13'  N 

176°28’W 

0.  77 

0.  05 

10°00'N 

170°00’E 

-2.  16 

2.  19 

180°00' 

0.  00 

2.  19 

170°00'W 

2.  16 

2.  19 

0°00' 

170°00'E 

-2.  19 

0.  00 

180°00' 

0.  00 

0.  00 

170°00'VV 

2.  19 

0.  00 

10°00'S 

170°00'E 

-2.  16 

-2.  19 

180°00' 

0.  00 

-2.  19 

170°00'W 

2.  16 

-2.  19 
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Figure  4.  Azimuthal  Equal-Area  Projection  Centered  on  Boston,  Massachusetts. 
The  map  extends  55.  (>  km  north-south  and  82.  2  km  east-west.  Each  cell  of  the 
square  grid  corresponds  to  243.4  km^  of  the  earth's  surface 


Figure  5.  Azimuthal  Equal-Area  Map  of  the  South  Pacific, 
from  10°S  to  10°N,  170°E  to  170°W,  Centered  on  the  Inter¬ 
section  of  the  Equator  and  the  International  Datpline 


'1 

I 
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Figure  (>.  Meridians  and  Parallels  for  Germany 


Figure  7.  Meridians  and  Parallels  for  Lcuador. 
Point  I,  is  at  Lasso,  point  Q  at  Queto 
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4.  EQl  Al.-AREA  MAPPING  CENTERED  ON  Till  NORTH  POLE 


4.1  Mapping  Nortli  al  Latitude  Op 

(Vhen  the  central  point  (b  ,  A  )  is  taken  at  the  North  Pole  the  projection  heroines 

o  •  » 

the  well-known  Polar  Azimuth  .1  Equal-Area  Projection  for  the  northern  hemisphere. 
As  in  Eq.  (4)  the  arc  (i  )  becoir.es  the  co-latitude,  so  tluit,  from  Kq.  (21  the  map 
distance  (r)  from  the  north  pole  n>  ;i  point  (P)  becomes: 


=  ^2(1  -  sin  b)  ■ 


(Rs> 


(10) 


Since,  at  the  north  pole,  A  does  not  have  a  unique  value  as  defined  so  far,  it  needs 
to  be  redefined  as  a  reference  meridian  from  which  others  are  measured.  The 
Cartesian  coordinates  become: 


x  =  r  ■  cos  (A  -  A  )  (11) 

o 

v  -  r  •  sin  (A  -  A  ) 

o 


4.2  Polar  Equal-Area  Mapping  South  of  Latitude  by 

Equations  (11)  are  used  in  the  construction  of  equal-area  maps,  suitable  out¬ 
ward  from  the  central  point  as  far  as  the  full  northern  hemisphere,  as  is  done  in 
the  DMA  map  (Figure  3).  For  mapping  beyond  the  hemisphere,  or  at  some  other 
desired  distance  from  the  central  point,  a  new  approach  is  adopted. 

The  construction,  as  described  below,  pertains  primarily  to  the  southern 
hemisphere  of  an  equal-area  map  centered  on  the  north  polo.  While  it  may  be 
convenient  to  begin  this  mapping  at  the  equator,  it  could  be  started  at  some  latitude 
(by)  north  of  the  equator.  Such  mapping  will  produce  a  top-of-the-world  view  of 
the  continents  (for  example  Figure  3  or  Figure  10  or  Figure  13).  Kxoept  for 
Antarctica  the  continents  need  not  be  split,  since  the  large  land  masses  of  the 
southern  hemisphere  can  fall  into  three  of  four  quadrants. 

The  plotting  is  done  by  Cartesian  (x,  y)  coordinates  (Figures  11  or  14)  with  the 
origin  (0)  representing  the  north  pole  and  the  X-  and  Y-axes  representing  two 
meridians  at  right -angles  to  each  other.  A  key  meridian  is  drawn  as  a  diagonal 
straight  line  <C)Q)  where  Q  represents  the  south  pole,  one  of  four  such  points,  each 
in  line  of  the  four  quadrants. 

If  the  angle 

(3  -  XOQ  (12) 


24 


id  chosen  creator  than  45°,  the  map  will  be  longer  than  it  is  wide.  If  .i  -  4a  .  then, 
as  in  Figure  3  or  10,  the  map  will  fit  neatly  into  a  square.  If  a  -  aa  .  as  in 
Figure  13,  the  map  will  fit  neatly  onto  an  8-1  3  X11  page. 

In  Figure  11  or  14,  the  meridian  (X)  that  makes  an  angle  (■> )  w  uh  the  kev 
meridian  (OQ)  is  drawn  straight  from  <)  to  I,  where  T  is  on  tile  circle  representing 
the  parallel  (o  ).  The  relation  between  X  and  -j  is 

1  .1  -  X  +  X  (13) 

where  X  ^  is  the  longitude  of  the  X-axis.  In  the  DMA  map  X  ^  -  -30°  for  the  meridian 

30''\V.  Beyond  T  the  meridian  is  represented  by  a  curve  (TQ),  so  that  the  area 

(( iTQ)  between  the  two  meridians  <(  )SQ,  OTQ)  on.  the  map  is  proportionally  equal  to 

the  corresponding  area  of  the  globe. 

■> 

1'he  total  global  area  is  4ttK~.  Between  two  meridians  subtending  angle  •% 

*) 

(radians)  at  either  pole  tlie  area  is  2->K“.  on  the  map  the  corresponding  area  is 
•> 

2>(Rs)“.  But  from  the  north  pole  southward  to  latitude  (©._)  the  global  area  is 
•>  1 
y(l  -  sin  •  H“.  On  the  map  the  cor  responding  area  is 

•> 

area  t'TM  y(l  -  sin  o •  (Ks)"  (14) 

For  the  balance  of  this  development  assume 

Us  1 

and  write 


area  OTX  "  l(l  -  sin  ©  .) 


(la) 


For  any  point  in  the  northern  hemisphere  with  spherical  coordinates  (o,  X)  where 
o  >  the  distance  (r)  from  the  north  pole,  on  the  map,  is 


r  2(1 


sin  o) 


ll'il 


The  area,  A^(y),  enclosed  by  latitude  (o)  subtended  by  the  apex  ancle  (•>)  i.- 


A  ,  ( i )  a  ( 1  -  sin©) 
o 


(17) 


The  map  distance,  (OT),  is  given  by 


1  T  J  2(1  *  s>n  °T) 

(18) 

Set 

CtQ  k 

(IP) 

The  straight  line  VT  is  drawn  tangent 
and  QV  is  drawn  perpendicular  to  VT, 

to  the  circle  representing  the  latitude  (p  ), 
making 

VT  =  k  sin  y 

(20) 

VQ  k  cos  •>  -  r 

(21) 

Area  \TQ  -  area  OTQ  -  area  OTN 

CO 

1 

1 

(f. 

5 

o- 

■  -y(l  +  sin  b-j.) 

(22) 

Also, 

Area  \TQ  Area  OTS  -  area  ()T\  +  area  VTQ  -  area  VSQ 
Hence, 

•y(l  +  sin  <>T)  =  1/2  r^,  tan  y  -  -yd  -  sin  $  )  +  A^  -  1/2  VQ2  tan  y 
where 

Ac  =  area  VTQ  (23) 

After  re-arrangement  of  terms, 

Ac  -  2-y  +  1/2  tan  7  •  (VQ2  -  r2  )  (24) 

At  this  stage  a  decision  must  be  made  on  the  shape  of  the  curve  TQ.  Handily,  if 
arbitrarily,  it  is  made  either  elliptical  or  hyperbolic.  The  two  alternatives  are 
described,  and  used  here. 
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The  frame  of  reference  (Figure  11)  is  transferred  f.  the  (X,  \  t-plnno  t.  ■  the 
(X1,  Y')-p)a»e  where  the  X'-axis  is  parallel  to,  but  not  identical  •«.  ith,  nr  irte  VQ, 
and  the  V'-axis  is  taken  along  the  tangential  line  VT.  l.et  the  distance 

O'T  =  b  (2'i) 

For  a  hyperbolic  shape  the  equation  for  the  curve  TQ  can  be  a  ritten 

x'  =  a  y'2  -  b“  +  g(y'  -  b)  (271 

(y  »•  b) 

where  the  parameters  a,  b,  g  need  to  be  determined  to  conserve  the  equal-area 
characteristic  of  the  map. 

For  Q,  (x1  -  VQ,  y  '  -  b  +  VT),  Fq.  (27)  gives 


VQ  a  /(b  +  VT)2  -  b2  t  g  .  VT 


Differentiating  Fq.  (27),  we  obtain 


1  r  (ay’/'^jy'2  -  b2  +  g)  •  (dy'/dx1) 


lienee,  at  T  (x'  =  0,  y'  -  b),  dy'/dx'  -  0 


Or,  the  slope  of  the  curve  TQ  makes  it  a  smooth  continuation  of  the  straight  line  OT 
into  the  southern  hemisphere,  a  desirable  effect. 
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At  Q  make 


angle  TQV  =  2y ' 

leaving  the  actual  value  of  y'  temporarily  undetermined.  At  Q 
dy'/dx'  =  tan  2-y1 

which,  together  with  Eq.  (29)  gives 
tan  2y(  =  h/(ac  +  gh) 
where 

c  -  (b  +  VT)/b 
h  -  ^/c2  -  1 

which  together  with  Eq.  (28)  give  expressions  for  a,  g  in  terms  of  b,  thus: 
a  =  h  ■  (VQ  -  VT  •  cos  2y')/VT 
g  =  cot  2-y'  -  ac/h 

The  area  VTQ  given  by  Eq.  (24)  is  also  given  by 
b+VT 

A'c  --  J  x'  dy' 

y=b 

which,  with  Eq.  (27)  becomes 

A'c  =  ab2[c  .  h  -  ln(c  +  h)i  /2  +  g  .  VT2/2 

So  far,  the  quantities  b  and  y'  have  been  assumed  known.  But  they  must  have 
that  will  provide  solutions,  such  that 


(SO) 


(31) 

(32) 

(33) 


(35) 

values 

(3  fi) 
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With  some  testing,  •>'  was  chosen  to  be  somewhat  larger  than  •>,  thus: 


+  ^/f2  -  (22.5  -  0>2  -  rj f2  -  (22.  5) 2 


where  f;  -  50  has  been  found  to  provide  viable  solutions.  The  value  of  b  is  found 
as  accurately  as  desired  by  trial  and  error,  to  satisfy  Eq.  (30).  With  b  known, 
the  values  o."  c,  h,  a  and  g  are  obtained  from  Eqs.  (32),  (33)  and  (34).  (See 
Table  4.  ) 


Table  4.  Table  of  Parameter  Values  for  the  Square  Equal-Area  Map 
Centered  on  the  North  Pole,  with  Ilyperbolically  Shaped  Meridians  in 
Southern  Hemisphere 


0 

VQ 

VT 

a 

b 

g 

0° 

0.3582 

1.  7725 

2. 1365  5 

1.  677  -6 

-2.  1365  5 

10° 

0.  03 91 

1.4377 

3. 6518 

0.  0934 

-3. 4374 

20° 

0.  857(1 

1.  0593 

4. 6118 

0.  0594 

-4.  0542 

30° 

1.  0070 

0.  6488 

6.  8421 

0.  0341 

-5.  6404 

O 

O 

1.  0829 

0.  2185 

5.  1095 

0.  0553 

-1.3144 

44. 99° 

1.0924 

0.4375  -4 

7.4164  2 

0.  7254  -4 

9.  5618  2 

For  each  y  the  curve  TQ  can  be  plotted  by  coordinates  (x',  y1)  using  Eq.  (27). 
Actually  it  is  best  to  find  the  (X,  Y)-coordinates  for  map  plotting: 

x  =  (x1  +  \^2)  cos  0  -  (y1  -  b)  sin  0  (3  8) 

y  =  (x'  +  -J2)  sin  0  +  (y'  -  b)  cos  0  (39) 

where 

0  =  P  -  y  (40) 

The  above  solution  yields  pairs  of  coordinates  (x,  y)  for  points  on  a  meridian 
(X),  But  this  has  left  the  problem  of  finding  (x,  y)  to  correspond  to  a  given  latitude 
(<t>)  as  well  as  longitude  (A). 
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Since  the  hyperbolic  shape  of  the  meridians  in  the  southern  hemisphere  is 
assigned,  the  shapes  of  the  parallels  must  allow  as  a  consequence,  subject  to  the 
constraint  that  the  map  is  equal -area.  Consider  the  point  1'  (Figure  12)  on  the 
hyperbolic  curve  TQ  representing  the  meridian  .if  longitude  (X).  The  point  F 
corresponds  to  the  earth  point  (b,  X)  where  d>  is  negative  in  the  southern  hemisphere. 
If  an  adjacent  curve  (T^Q)  is  drawn  such  that  the  angle  subtended  by  T^T  at  (l  is 
small  (67),  then  the  small  incremental  area  (F^FQ)  in  the  southern  hemisphere 
must  be  made  equal  to  the  corresponding  small  increment  in  the  northern  hemi¬ 
sphere.  That  is, 

area  FjFQ  =  67(1  -  si (41) 

But 

n 

area  FjFQ  -  f  (x  -  x^)  •  dy  J  £  (x  -  Xj).  •  6y  (42) 

y-y W  i=l 


where  6v  is  a  finite,  but  ronstant,  increment  that  we  ran  make  as  small  as  we 
please  by  dividing  the  Y-distance  from  Q  to  T  into  as  many  equal  parts  as  we 
please.  We  make 

by  --  U'tr  -  s2  •  sin  9)/ N  (43) 


where  N  is  large  (1000  say).  The  summation  of  terms  in  Kq.  (42)  is  such  that  the 
sum  from  Q  to  F  is  equal  to  the  area  given  by  Eq.  (41).  The  number  of  such  terms 
(n)  increases  as  the  latitude  decreases  in  magnitude. 

For  the  i1*1  term  in  Eq.  (42) 

v  -  -  i*y  (44) 

With  x1  given  in  terms  of  y1  in  Eq.  (27),  then  Eq.  (39)  yields  a  value  for  y'  as 
follows: 


y’ 


(b  +  s) 


(45) 


where 


P  =  (y  -  s/2  sin  9) fa.  sin  6 


(46) 
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The  value  of  x  is  obtained  on  the  curve  QFT.  Likewise  the  value  of  Xj  is  ob¬ 
tained  on  the  curve  QF^T^.  The  summation  of  successive  terms 

(x  -  x,)^  •  6y  o  <  i  <  n 

gives  the  approximate  integration  in  Eq.  (42)  for  the  area  FF^Q.  The  summation 
increases  as  the  magnitude  of  <j>  decreases. 

An  alternative  development  was  used  when  9  was  very  small  and  for  latitude 
close  to  the  equator.  Instead  of  the  area  F^FQ  the  complementary  area  was  used: 

Area  T.TFF^  =  iy  •  sin  (- <j> ) 

for  which,  also: 

x(<£)  n 

Area  T ^TFFj  =  J  (y1  -  y)  •  dx  =  ^  (y}  -  y).  •  6x  (49) 

x=>/2  cos  8  1  * 


where 


6x 

=  («/tr  -  s' 2  •  cos 

0)/N 

(50) 

For  the 

.th 

i  term 

x  = 

sr2  c  os  $  +  i  ■  ft 

X 

(51) 

With  x  thus  given  and  x' 

given  in  tei 

rms  of  y<  in  Eq. 

(27)  the  value  of  v'  is  found  from 

y'  ; 

-■  (b  +  S)  - 

r  2,,  2 

-  (P  /<P  - 

1)} 

(52) 

where 

P  = 

(sin  0  -  g  •  cos 

0)/( a  •  cos 

8) 

(53) 

P  ; 

(x  -  sr2  •  cos  8) /(a  •  cos  9) 

(54) 

s 

(b  -  p  •  p)/(p2  - 

■  1) 

(55) 

With  both  x',  y'  thus  found,  the  value  for  y  is  obtained  from  Eq.  (39). 
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Likewise  the  value  of  y1  is  obtained.  The  summation  of  successive  terms 

(y j  -  y)  •  fix  o  <  i  <  n 

gives  the  approximate  integration  in  Eq.  (49)  for  the  area  T^TFF^.  The  summation 
increases  as  the  magnitude  of  ^  increases. 

The  above  procedure  gives  points  in  the  lower  half  of  the  upper  right-hand 
quadrant  of  the  southern  hemisphere,  below  line  OQ  (Figure  11  or  12).  For  the 
map  above  OQ  the  (X,  Y)-coordinates  are  as  follows: 

For  45  <  (X  -  X  )  £  90°  x($,  X)  =  y(d>,  90  -  X  +  Xq)  (56) 

y($,  X)  =  x($,  90  -  X  +  Xq) 

For  points  in  the  other  three  quadrants: 

For  90  <  (X  -  X  )  <  180  x(<t>,  X)  =  180  -  X  +  >  )  (57) 

o  o 

X)  -  y (<t>,  180  -  X  +  x  ) 

For  -180  <  (X  -  Xq)  <  -90  x(</>,  X)  =  -x(<t>.  180  +  X  -  Xq) 

y (<t>,  X)  =  -y (4>,  180  +  X  -  Xq) 

For  -90  <  (X  -  Xq)  <  0  x(<j i,  X)  =  x(0,  -X  +  X^) 

y(<t>,  X)  =  -y (<t>,  -X  +  X^) 

Since  the  map  is  equal-area,  it  is  possible  to  subdivide  it  into  cells  of  equal 
area,  that  correspond  to  likewise  equal  areas  on  the  globe,  by  latitudinal  and 
longitudinal  lines  or  curves.  But,  additionally  since  the  map  is  square,  it  can  be 
divided  into  equal  areas  by  a  square  grid  superimposed  on  the  map. 

For  computer  programming,  the  above  procedure  should  be  followed  in  steps 
(Appendix  B).  The  Cartesian  coordinates  (x,  y)  corresponding  to  global  points  at 
<t>  -  0(-10°)-90°,  X  =  0(10)40,  44.  99  are  shown  in  Table  5. 

4.  2.  2  ELLIPTIC  ALLY  SHAPED  MERIDIANS  SOUTH  OF 
LATITUDE  dT  (Figure  13) 

In  Figure  14  the  frame  of  reference  is  transferred  from  the  (X,  Y)-axes  to  the 
(X1,  Y')-axes.  The  X'-axis  is  taken  along  the  line  VT.  The  Y'-axis  is  parallel  to, 
but  not  identical  with,  the  line  VQ.  Let  the  distance 

O'T  =  a  (58) 
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Table  5.  The  (x,  y)-Coordinates,  Corresponding  to  Points  (i j>,  A.)  on  the  Global 
Surface,  on  a  Square  Equal-Area  Map,  Which  is  Lambert  Azimuthal  in  the 
Northern  Hemisphere  and  With  Hyperbolic  Meridians  in  the  Southern  Hemi¬ 
sphere  (Figure  10).  Assume  X  =  0 


l.at 

-  4,  0 

-10 

-20 

-30 

-40 

-50 

-fiO 

-70 

-80 

-90 

Long 

Kquator 

io“s 

20°S 

30°S 

South 

Pole 

X 

0° 

X 

1.414 

1.  772 

1.  772 

1. 772 

1.  772 

1.  772 

1.  772 

1.  772 

1.  772 

1. 772 

y 

0.  000 

- 

- 

- 

- 

- 

- 

- 

- 

1.  772 

10 

1. 393 

1.  50  ti 

1 . 500 

1.  (Ml 

1.  (598 

1.  722 

1.  738 

1.  751 

1.  703 

1.  772 

0.  24 ‘i 

0.  270 

0.  334 

0.444 

0.  003 

0.  804 

1.027 

1.  2(>9 

1.  519 

1.  772 

20 

1.329 

1. 43  5 

1.  o20 

1.  570 

1.  017 

1. 052 

1.  084 

1.  714 

1. 743 

1.  772 

0.  484 

0.  535 

0.  020 

0.  744 

0.  894 

1.058 

1. 231 

1. 409 

1.  589 

1.  772 

30 

1.  225 

1.322 

1.397 

1.458 

1.  512 

1. 505 

1.  010 

1.  Mi  8 

1. 720 

1.  772 

0.  707 

0.  775 

0.  801 

0.  907 

1.  080 

1.  210 

1.  351 

1. 489 

1.  030 

1.  772 

40 

1.083 

1.  171 

1. 249 

1.323 

1.397 

1.471 

1.  545 

1.  020 

1.  090 

1.  772 

0.  90!) 

0.  990 

1.  073 

1.  105 

1. 258 

1.350 

1. 450 

1.  500 

I.  005 

1.  772 

44. 

99 

X 

1.  000 

1.085 

1.  109 

1.253 

1.338 

1. 425 

1. 511 

1.  509 

1. 080 

1.  772 

y  - 

1.000 

1.085 

1.  108 

1. 253 

1.  338 

1.  425 

1.  511 

1.  599 

1.  080 

1.  772 

The  parameters  (1,  k,  0,^,  and  consequently  r.p,  have  presently  unassigned  values. 
For  an  elliptical  shape  the  equation  for  the  curve  IQ  can  be  written: 

y‘/b  =  *  <x'/a)2  -  g  •  a  .  (1  -  x'/a)  (59) 

where  the  parameters  a,  b,  g  need  to  be  determined  to  conserve  the  equal  area 
characteristic  of  the  map. 

For  the  point  Q  (x1  =  a  -  VT,  y'  -  VQ)  Eq.  (59)  gives 

VQ  /  b  /syi  -  ((a  -  VT)  a  -  g  •  VT  (HQ) 


Differentiating  Eq.  (59),  we  obtain 
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Figure  14.  To  Illustrate  the  Construc¬ 
tion  of  the  World  Map  With  Elliptieallv 
Shaped  Meridians  in  the  Southern  Hemi¬ 
sphere 


Hence,  dy'/dx1  —  as  x'  -»  a  or,  the  slope  of  the  curve  TQ  at  T  is  perpendicular 
to  VT,  in  order  to  make  the  curve  TQ  a  smooth  continuation  of  the  line  OT,  a 
desirable  effect. 

At  point  Q,  make 

<  TQV  =  2y'  (02) 

where  -y'  is  slightly  larger  than  y.  Arbitrarily  we  make 

V  =  7  •  exp  {y1/32  •  (3  -  y)/3  }  (63) 

where  the  angular  measurements  in  the  bracket  are  in  radians,  which  makes 
y'  =  y  at  y  =0 
>  y  at  0  <  y  <  P 
=  7  at  y  =  3  . 
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dy'  /  dx'  -  -  cot  (2-> ') 


which,  together  with  Eq.  (fil),  gives 


-cot  (2y')/b  =  g  -  {1  -  (VT/a)l/aJl  -  {(a  -  VT)/a|‘ 


Equations  (tiO)  and  (t>4)  can  be  solved  to  obtain  expressions  for  b,  g  in  terms  of  a. 


b  =  h  •  a  [  VQ  -  VT  ■  cot  (2y )]  /VT 


g  =  h/VT  -  VQ/b  •  VT 


where 


h  =  Jl  -  {(a  -  VTl/af 


The  area  VTQ  given  bv  Eq.  (24)  is  also  given  by 


V  •  dx' 


x'  =  a-VT 


which,  with  Eq.  (59),  becomes 

A'  =  (ab/2)[tr/2  -  (1  -  VT/a)  •  h  -  sin'1  (1  -  VT/a)  -  g  •  a  (VT/a)2)]  (09) 

Solution  for  a  is  obtained  by  making 

A’.  =  Ac  (70) 

This  can  be  done  by  trial  and  error,  as  accurately  as  desired.  It  has  been  done  to 
12  decimal  places  in  approximately  40  iterations.  When  the  value  of  a  is  put  in 
final  form,  the  values  of  b  and  g  are  given  by  Eqs.  (05)  and  (tifi). 

Once  the  values  for  a,  b,  g  are  obtained,  for  the  given  •>,  the  curves  TQ  can 
be  plotted  by  coordinates  (x',  y')  through  Eq.  (59).  Actually  it  is  best  to  find  tile 
(x,  y)-coordinates  for  map  plotting: 
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x  =  (v1  +  r^)  sin  8  +  (x1  -  a)  cos  6 


y  r  (y1  +  r^,)  cos  0  -  (x<  -  a)  sin  8 


0  -  n  2  -  ,i  +  -y  (7:0 

The  foregoing  solution  yields  pairs  of  Cartesian  coordinates  for  plotting  merid¬ 
ians.  Since  the  elliptical  shape  of  the  meridians  is  assigned,  the  shapes  of  the 
parallels  must  follow  as  a  consequence,  subject  to  the  constraint  that  the  map  is 
equal-area. 

Consider  the  point  F  (Figure  15)  on  the  elliptical  curve  TQ  representing  the 
meridian  of  longitude  (A.).  The  point  F  corresponds  to  the  earth  point  (b.  A.)  where 
d  is  negative  in  the  southern  hemisphere.  If  an  adjacent  curve  (TjQ)  is  drawn  such 
that  the  selected  angle  subtended  by  T^T  at  o  is  small  (6y),  then  the  small  incre¬ 
mental  area  (FjFQ)  in  the  southern  hemisphere  must  be  made  equal  to  the  corre¬ 
sponding  small  increment  in  the  northern  hemisphere  as  given  by  Eq.  (17).  That 


area  F.FQ  =  {1  -  sin  (-d)f 


k  sin  .1 

FiFQ  ,  / 

v-y(b) 


(x  -  Xj)  •  dv 


2 (x  -  xi) 


where  6y  is  a  finite,  but  constant,  increment,  which  we  can  make  as  small  as  we 
please  by  dividing  the  v-distance  from  Q  to  T  into  as  many  equal  parts  as  we  please; 
that  is,  we  make 

6 y  f(k  sin  ,1  -  r.j.  sin  (,i  -  y)};\  (7*0 

where  N  can  be  made  large,  1000  say.  The  sum  nr  lion  of  terms  in  Fti.  (75)  is  tr  dc 
until  th“  sum.  from  Q  to  1  ,  is  equal  to  the  area  given  by  Kq.  (74).  The  number  of 
such  terms  ( n >  increases  is  the  latitude  to)  decreases  m  magnitude.  I  lie  co.e  !- 
mates  (x,  y)  at  which  Kq.  (75)  is  satisfied  become  the  loc  ma  s  of  the  point  that  has 
sfihc  i  ■ ; j i  i  |.,,r  Ima'cs  ( ■  0,  A). 


Q 

/I 


Figure  15.  To  Illustrate  the 
Determination  of  the  (X.5:)- 
coordinates  on  the  Elliptical 
Curve  (TQ)  to  Correspond  to 
Spherical  Coordinates  (<Ji,  X) 


(77) 


(78) 


and  where 


q  -  ■  (y  +  m) 

With  y1  given  by  Eq.  (59),  then 

x  =  (y1  +  r^)  sin  6  +  (x'  -  a)  cos  0  (79) 

This  value  of  x  is  obtained  on  the  curve  QFT.  Likewise  the  value  of  is 
obtained  on  the  curve  QFjT^. 

The  summation  of  successive  terms 

(x  -  x,)^  •  6y 

gives  the  approximate  integration  for  area  F^FQ.  The  summation  increase's  as  the 
magnitude  of  <j>  decreases.  Sufficient  accuracy  is  obtained  when  is  sufficiently 
small,  although  it  must  not  be  too  small  for  computer  operability. 

The  above  procedure  was  used  to  produce  the  coordinates  in  the  upper  right  - 
hand  quadrant  of  the  map  (Figure  13),  at  the  intersection  of  meridians  and  parallels 
at  10°  intervals  (Table  0).  The  steps,  for  the  computer  are  given  in  Appendix 

Several  key  points  and  additions  to  the  above  analysis  are  needed.  For  the 
elongated  map  (3  (Figure  14)  is  greater  than  45°.  The  lower  part  of  the  upper  right- 
hand  quadrant,  below  line  OQ,  is  derivable  by  direct  application  of  the  above  pro¬ 
cedure.  For  the  upper  part,  above  line  (IQ,  the  computations  should  be  done  with 
the  mirror  image,  in  which  ;3  is  less  than  4.5°  and  giv  -n  by 

|3m  =  90°  -  ,3  (80) 

If  the  Cartesian  coordinates  so  obtained  are  designated  (x  ,  v  ),  then  the  primary 
coordinates  are  given  by  the  equations: 

x  =  x^  cos  2(3  +  ym  sin  2(3  (81) 

y  =  ym  sin  2(3  -  ym  cos  2(3 

Another  key  point  concerns  the  latitude  (^^,)  at  which  each  meridian,  beginning 
at  the  north  pole,  ceases  to  be  represented  by  the  straight  line  and  begins  its  ellip¬ 
tical  shape  ending  at  the  south  pole.  In  the  example  (Figure  13)  ^(X)  is  at  50°N 
along  the  X  =  axis  and  i^rp(Y)  at  20°N  along  the  Y-axis,  and  is  a  function  of  the 
longitude  in  the  remainder  of  the  quadrant,  as  follows: 
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Table  (i.  The  Rectangular  Coordinates  (x,  v)  of  the  Intersection  of  Meridians  and 
Parallels  at  10  Intervals  in  the  First  Quadrant  of  the  Klongated  Polar  Kqual -A rea 
Map 


.ongitude 

> 

X 

Latitude 

20°\V 

55° 

-20 

x  y 

10’ 

45 

-10 

X 

\v 

y 

3 

X 

0° 

5 

0 

y 

10l 

25 

10 

X 

y 

20°  i: 

15 

20 

x  y 

N  -  Pole 

0.  000 

0.  000 

0.  000 

0.  000 

0.000 

0.  000 

0.  000 

0.  000 

0.  000 

0.  000 

80°  N 

0.  174 

0.  000 

0.  172 

0.  030 

0.  1 04 

0.  080 

0.  151 

0.  087 

0.  134 

0.  112 

70 

0.347 

0.  000 

0.  342 

o.  otto 

0.326 

0.  119 

0.  301 

0.  174 

0.  266 

0.  223 

co 

0.  518 

0.  000 

0.  510 

0.090 

0.  481) 

0.  177 

0.  448 

0.  259 

0.  397 

0.  333 

50 

0.  084 

0.  000 

0.  074 

0.  119 

0.  043 

0.  234 

0.  592 

0.342 

0.  524 

0.440 

40 

0.  850 

0.  014 

0.  862 

0.  161 

0.820 

0.  302 

0.  732 

0.  42.3 

0.  647 

0.  543 

30 

1. 028 

0.  056 

1.  009 

0.  220 

0.  956 

0.  375 

0.879 

0.  523 

0.  777 

0.  659 

20 

1.  178 

0.  123 

1.  149 

0.  3  03 

1. 084 

0.  409 

0.  998 

0.  627 

0.  881 

0.  770 

10°  N 

1.320 

0.  218 

1. 279 

0.  411 

1. 203 

0.  384 

1.  103 

0.  74  5 

0.  979 

0.  890 

Equator 

1. 447 

0.  338 

1.  393 

0.  538 

1.309 

0.  715 

1. 200 

0.  873 

1. 070 

1. 017 

10°S 

1. 557 

0.484 

1. 493 

0.  080 

1.401 

0.  880 

l.  288 

1.  012 

1.  1  .4 

1.149 

20 

1.  848 

0.  853 

1.  577 

0.  855 

1. 479 

1.018 

l .  3  83 

1 .  1 59 

1. 233 

1.  287 

30 

1.  710 

0.  843 

1.  641 

1.0.3  7 

1 . 543 

1.  188 

1.431 

1.315 

1.  308 

1. 430 

40 

1.  787 

1. 054 

1.  687 

1.  2.31 

1 .  592 

1 .  3  83 

1.487 

1.  474 

1.375 

1.  374 

50 

t  791 

1. 284 

1.  713 

1.437 

1.  627 

1.  548 

1.  53  5 

1.  639 

1. 43  7 

1.  720 

80 

1.  790 

1.  527 

1.  720 

1. 652 

1.  647 

1.  738 

1.  571 

1.  805 

1.  493 

1. 870 

70 

1. 782 

1.  732 

1.  707 

1.  870 

1. 653 

1.  930 

1.  598 

1. 976 

1.  542 

2.  018 

80°S 

1. 707 

2.  046 

1.  674 

2.094 

1,  644 

2.  123 

1.815 

2.  146 

1. 586 

1.  167 

S-Pole 

1.622 

2.  317 

1.  622 

2.317 

1. 622 

2.  317 

1.  822 

2.  317 

1.  622 

2.317 
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Table  li.  The  Rectangular  Coordinates  (x,  y)  of  the  Intersection  of  Meridians  and 
Parallels  at  10  Intervals  in  the  First  Quadrant  of  the  Klongated  Polar  Kqual-Area 
Map  (Continued) 


Longitude 

> 

X 

Latitude 

3o"i: 

5 

30 

x  y 

40l 

-5 

40 

X 

*K 

y 

50°  i: 

-15 

50 

x  y 

(50 

-25 

00 

X 

'i: 

V 

70 

-35 

70 

X 

L 

V 

N-Pole 

0.000 

0.000 

0.  000 

0.  000 

0.000 

0.  000 

0.  000 

0.  000 

0.  000 

0.  000 

80".\ 

0.  112 

0.  134 

0.  087 

0.  151 

0.  080 

0.  104 

0.030 

0.  172 

0.  000 

0.  174 

70 

0.  223 

0.  2t>(» 

0.  174 

0.  301 

0.  119 

0.  320 

0.  000 

0.  342 

0.  000 

0.  347 

tiO 

0.  333 

0.  397 

0.  259 

0.  448 

0.  177 

0.  48(5 

0.  090 

0.  510 

0.  000 

0.  518 

50 

0.440 

0.  524 

0.  342 

0.  592 

0.  234 

0.  (543 

0.  119 

0.  074 

0.  000 

0.  084 

40 

0.  543 

0.  847 

0.  423 

0.73  2 

0.  289 

0.  794 

0.  147 

0.  832 

0.  000 

0.  84  5 

30 

0.  053 

0.  780 

0.  509 

0.  881 

0.342 

0.  940 

0.  174 

0.  984 

0.  000 

1.  000 

20 

0.  744 

0.  897 

0.  587 

1. 008 

0.  407 

1.090 

0.  204 

1.  132 

0.  000 

1.  147 

10°N 

0.834 

1. 020 

0.  870 

1.  129 

0.485 

1.  218 

0.  275 

1.  271 

0.  054 

1.  304 

Kquator 

0.  921 

1.  144 

0.  757 

1.253 

0.  575 

1. 343 

0.  370 

1. 401 

0.  158 

1.  443 

to°s 

1.008 

1.  272 

0.  848 

1.377 

0.  li 74 

1.485 

0.  482 

1.  527 

0.  284 

1.  574 

20 

1.093 

1.  403 

0.  941 

1.  502 

0.  780 

1.  580 

0.  (505 

1.  849 

0.  428 

1.  899 

30 

1.  178 

1.  535 

1.037 

1.  028 

0.  891 

1. 704 

0.  73  5 

1.  785 

0.  580 

1.  818 

40 

1.  258 

1.  087 

1.  135 

1.  749 

1.008 

1. 818 

0.  873 

1. 878 

0.  742 

1.925 

50 

1. 337 

1.  800 

1.  233 

1.  870 

1.  125 

1. 929 

1. 017 

1.  981 

0.  910 

2.  026 

00 

1.414 

1,  932 

1.332 

1.988 

1.  247 

2.03  8 

1.  183 

2.  079 

1.  083 

2.  118 

70 

1.487 

2.093 

1.  430 

2.  103 

1. 370 

2.  137 

1.312 

2.  109 

1. 239 

2.  199 

80°S 

1.  557 

2.  191 

1.  527 

2.213 

1.  539 

2.  018 

1.405 

2.  249 

1.438 

2.  287 

S-  Pole 

1.  822 

2.  317 

1.  922 

2.317 

1.  822 

2.317 

1.  022 

2.  3 1 7 

1.  022 

2.  .3  1 7 

d>T  =  50  -  1/3  (X  -  \q) 

More  generally 

<*>T  =  <^T(X)  -  (X  -  Xq)  .  (^T(X)  -  0t(Y»/9O  (82) 

For  the  other  three  quadrants  of  the  map,  the  (x,  y) -coordinates  can  be  found 
in  terms  of  corresponding  positions  in  the  first  quadrant  by  the  transformations  in 
Eq.  (57), 

Since  the  map  is  equal-area,  it  is  possible  to  subdivide  it  into  cells  of  equal 
area  that  correspond  to  equal  areas  on  the  globe  as  well  (Figure  16).  For  the 
division  of  each  hemisphere  into  M  equal  cells  by  latitudinal  and  longitudinal  di¬ 
viders,  the  latitude  down  to  the  nth  ring  will  be  given  by 

<l>  =  sin  *  (1  -  m  /M)  (83) 

’n  n' 

where  m_  is  the  number  of  cells  down  to  latitude  <t>  .  It  is  best  to  choose  M 
n  n 

divisible  by  four,  to  preserve  the  symmetry  of  the  four  quadrants. 

There  is  a  restriction  in  the  use  of  the  above  procedure  to  construct  an  equal- 
area  map.  The  solution  is  not  always  possible  with  the  arbitrarily  selected  values. 
Figure  13  was  obtained  after  several  tests  with  varying  values  for  the  parameters: 
tf-p,  k,  y',  before  the  values  were  selected  as  /i  =  55°,  k  =  2<.r2,  <J.j.(u)  -  aO°N, 
</>T(i)  ~  20°N  and  y1  and  <j>T  given  by  Eqs.  (63)  and  (82). 

The  Cartesian  coordinates  (x,  y)  corresponding  to  global  points  at 
<t>  =  90°(-10)-90°,  X  =  0°(  10)90°  are  shown  in  Table  6. 

5.  ILLUSTRATIONS 

Maps  of  mesoscale  characteristics  are  "everywhere,”  and  world  maps  of 
specifics  are  commonplace.  This  report  was  written  to  advance  the  use  of  equal  - 
area  maps,  both  large-scale  for  mesoscale  studies  and  small-scale  for  world-wide 
depiction.  Admittedly,  there  have  been  few  applications,  so  far,  with  which  this 
report  could  be  illustrated.  Consequently  there  may  be  problems  as  yet  unantici¬ 
pated. 

The  world  map  issued  by  the  Defense  Mapping  Agency  (PEA -3)  has  been  used 
satisfyingly.  Figure  17  shows  isopleths  of  average  cloud  cover,  throughout  the 
world,  in  January  at  local  noontime.  The  centers  of  high  frequency  of  cloud  cover 
are  marked  with  the  plus  sign  (t),  usually  better  than  8.' 10,  and  the  centers  of 
low  frequency  are  marked  with  the  minus  sign  (-),  often  lower  than  2/  10.  Figure  18 
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Figure  18.  The  Square  Equal-Area  Map  Showing  Average  Cloud  Cover,  in 
January,  Local  Noontime 


shows  the  same  information  on  the  square  equal-area  map,  on  which  there  actuallv 
are  no  breaks  in  the  isopleths  from  one  quadrant  to  the  next,  but  at  a  cost.  They 
are  crowded  in  the  higher  latitudes  of  the  southern  hemisphere. 
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6.  CONCLUDING  KKMARK 


Much  of  this  report  has  been  devoted  to  the  eons!  met  ion  of  tin  poko  qu  d-.o 
maps  of  the  world  inasmuch  as  they  are  new  "to  the  market."  However,  ’he  mam 
objective  is  to  advance  the  use  of  azimuthal  equal-area  mapping,  which  car.  be 
computerized  and  expedited  in  all  details. 
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Appendix  A 

Program  Steps  for  Plotting  Cartesian  Coordinates  (x,y)  Corresponding 
to  Geographic  Coordinates  (d> , A.)  on  an  Azimuthal  Equal-Area 
Projection  Centered  on  (<40.A0)  Given  the  Map  Scale  (s) 


Step  1.  Enter  the  radius  of  the  earth  :  R  -  2.  508  X  108  inches 
Enter  the  scale  (s) 

Step  2.  Enter  the  coordinates  of  the  central  point  (<j>  A  ) 

-90°  <  <4  <90°  -180°  <  A  <  180° 

o  o 

Note:  When  (4^  =  90°,  the  map  is  the  Polar  Azimuthal  Equal-Area  Map, 
which  is  treated  separately. 

Step  3.  Enter  the  coordinates  of  the  station  to  be  plotted  (<4,  A). 

Step  4.  Find  cos  v  -  sin  <4  sin  <4  +  cos  6  cos  <4  cos  (A  -  A  ) 

o  o  o 

Find  v,  considered  less  than  90°  and  positive. 

Step  5.  Find  sin  A  =  {cos  <4  •  sin  (A  -  A^jJ/sin  v 

Find  A,  which  will  be  negative  when  (A  -  A  )  is  negative. 

Note:  Choose  |a|  less  than  90°  at  this  stage. 

Step  fi.  Find  cos  <j>s  =  1  -  sin^  4  )\  { 1  -  sin^  <j>  •  sin^  (A  -  A^>  ) 
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Step  7. 
Step  8. 
Step  9. 


i 


Find  which  will  be  negative  when  is  negative,  as  it  is  in  the 
southern  hemisphere. 

If  4  <  $s,  replace  A  by  (180°  -  A) 


Find  r 


-  cos  v)  •  (Rs) 


Find  x  =  r  sin  A 
y  =  r  cos  A 

Print  and/cr  plot  (x,  y). 


i 
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Appendix  B 

Program  Steps  for  Plotting  Cartesian  Coordinates  (x,y)  Corresponding 
to  Geographical  Coordinates  (c> ,  X)  on  the  Square  Equal-Area  Map 
(Limited  to  the  Octant:  0  <  x  <  *,  0  <  y  «  x) 


Assume  that  the  earth's  radius  (I?)  times  the  map  scale  (s)  are  such  that 
Rs  =  1.  Or,  assume  that  all  map  distances  are  multiplied  by  (1  Rs). 

Step  1:  Enter  X  ,  the  meridian  of  reference  (-20°  in  the  DMA -map) 

r  o'  r 

Enter  ;J  =  45°  (for  a  square  map) 

Enter  k  =  n'2t 

Enter  <>X  =  0.001°  (a  suitable  value) 


Step  2. 
Step  3. 


Enter  ^  =  0  (at  the  equator) 

r.j,  v  2  (at  the  equator) 

N  -  1000  (or  other  suitable  large  number) 
?  -  50  (for  a  viable  solution  For  y') 

Enter  latitude  (<j>)  and  longitude  (X)  of  the  station. 
If  <j>  <  0,  go  to  Step  4 


For  ij>  0  (northern  hemisphere) 


find  r  -  / 2(1  -  sin  tj>) 


find  x  =  r  cos  (X  -  X  ) 

o 

find  y  -  r  sin  (X  -  X  ) 
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(END) 


Stop  4. 


Stop 
Stop  H. 

Step  7, 

Step  8. 

Step  9. 


Step  10. 


Stop  11. 


Fur  0*0 


1  f  0 


-l go  to  stop  1 1. 


For  6  <  -  1  :>l> 


find  Ao (6X)  6X{1  -  sin  (-o)l 

Enter  t  (x  -  x^  by  -  0  (initial) 
y  (initial) 

Find  Sv 


|N  tr  -  r  ■  sin  (X  -  A.  ) }  ■'  \ 


Use  SUBROUTINE  A  to  find  (and  save)  0,  b.  a,  g,  p. 

Find  X^  X  ■*-  fiX 

Use  SUBROUTINE  A  to  find  (and  save)  f)  b  a^,  g^,  p^ 
With  6,  b,  a,  g,  p,  y 

Use  SUBROUTINE  B  to  find  (and  save)  x 
With  0j,  b1,  a1.  pv  y 
Use  SUBROUTINE  B  to  find  (and  save)  Xj 
Find  (x  -  Xj)  •  6v 


and  find  new  ^  (x  -  x^).  •  by 


i  1 


If 


E(x  -  x^jfiy  <  A^(dy)  go  to  Step  10. 


For  E(x  -  x , ) .  6 v  >  A  (6v) 

I  I  o 

read  (x,  y)  (END) 

Decrease  y  by  6y  to  give 
y(new)  =  y(old)  -  by 
Go  to  Step  7. 


For 


-15 


find  A  ,(6X)  =  6X(sin  ( -<^)) 

o 

Enter  E(y^  -  v )6x  -  0  (initial) 

x  *  vr2  •  cos  (X  -  X  )  (initial) 


Step  13. 


X  •  6\ 


Kind  X 

Use  SUBROUTINE  A  to  find  (and  save)  0j,  bj,  a^,  g ^ ,  p  1 
Step  14.  With  f>,  b,  a,  g,  p,  x 

l  sc  Sl'RIU  it'TTNK  B  to  find  (and  save)  y 
Step  la.  With  Oy  bj,  a  gy  py  x 

Kse  SUBROUTINE  B  to  find  (and  save)  Vj 
Step  lti.  Kind  (Vj  -  y)  •  6x 

and  find  new  -  y)j  •  6x 

If  £(1^  -  y).  '  <5x  <  A^tty)  go  to  Step  17. 
if  t(yt  -  y), '  «x  S  A^&y) 
read  (x,  y)  (END) 

Step  17.  Increase  x  by  6x  to  give 
x(new)  =  x(old)  +  6x 
Go  to  Step  14. 

SUBROUTINE  A. 

(Using  previous  entries  including  \  find  0,  b,  a,  g,  p.) 
Step  1.  Find  0  -  X  -  X^ 

7  -  3  "  0 
VT  =  k  sin  y 
VQ  =  k  cos  y  -  i  y 

Ac  -  2y  (rad)  +  1/2  tan  y  (VQ2  -  r^.) 

T'°  =  7°  +  /s/?2  -  (22.5  -  e)2  -  ^/?2  -  (22.  5)2 
Step  2.  Initialize  b(u)  =  100 
bit)  =  0 

v  -  0 

Step  3.  Increase  p  by  1. 

Find  b  =  ^  b(u)  +  b(< ) }  /2 
c  =  1  +  VT/b 


a  =  h[VQ  -  VT  •  cot  2y<]  /VT 

g  -  cot  2y'  -  a  •  c/h 

AJ,  =  a  .  b2  [ch  -  l  n(c  +  h)J  /2  +  gVT2/2 

Step  4.  If  A'  >A 

c  c 

Make  b(i )  =  b 
Go  to  Step  5 

If  A'  <  A 

c  c 

Make  b(u)  =  b 
Step  5.  If  i/  >  38 

Go  to  Step  6 

For  v  <  38 
Go  to  Step  3 
Step  6.  If  $  <  -15° 

Find  p  =  (g  +  cot  6)/ a 

(RETURN) 

If  <t>  >  -15° 

find  p  =  (tan  6  -  g)/a 

(RETURN) 


SUBROUTINE  B. 

(Using  0,  b,  a,  g,  p,  y  or  x,  find  x  or  y.  ) 
Step  1.  If  <»  5  -15° 
go  to  Step  2 
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Step  2. 


X  -  (X1  +  r  >  cos  f)  -  (y  '  -  b)  sin  0 

(RETURN) 

For  >  -15° 

p  --  (x/ros  6  -  »•  |/a 

s  =  (b  -  p  •  p)/(p2  -  1) 

y'  =  (b  +  s)  +^/s2  +  {p2/(l  -  p2)} 

x'  -  a ^Jy'2  -  b2  +  g(y'  -  b) 

v  =  (x'  +  r^,)  sin  <j>  +  (y1  -  b)  ros  ij> 

(RETURN) 
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Appendix  C 

Program  Steps  for  Plotting  Cartesian  Coordinates  (x,y)  Corresponding 
to  Geographical  Coordinates  (4>  ,X )  on  a  Polar  Equal-Area  Map 
with  Elliptical  Meridians  South  of  Latitude  <t> T  (Limited  to  the 
Segment  0<x<«,  O^ysx  tan  3) 


O 

Assume  that  the  earth's  radius  (R  3  2.  508  X  10  inches)  multiplied  by  the  map 
scale  (s)  gives  Rs  =  1.  Or,  assume  that  all  map  distances  are  multiplied  by  (1/RS). 

Step  1.  In  this  step  the  bracketed  figures  apply  to  Figure  13. 

Enter  X  ,  meridian  of  reference  (-20°) 
o’ 

3,  azimuthal  angle  of  the  key  meridian  (55°) 

k,  length  of  the  key  meridian,  the  diagonal  straight  line  from 
north  pole  to  south  pole  (2  2 ) 

6X,  a  very  small  increment  cf  longitude  (0.001°) 

$,j,(X),  the  latitude  at  which  the  map  is  split  along  the  X-axis 
(50°N) 

4>^,(Y),  the  latitude  at  which  the  map  is  split  along  the  Y-axis 
(20°N) 

N,  a  large  number  (1000) 

Step  2.  Enter  latitude  (<)>)  and  longitude  (X)  of  the  station. 


til 


Step  3. 


Step  4. 


Step  5. 


Step  6. 

Step  7. 


Step  8. 


Step  9. 


Find  0T  =  «*T(X)  -  (X  -  Xq)  .  (*T<X)  -  4>t(Y))/90 
Find  rT  =  ^2(1  -  sin  $T> 

If  ij>  <  go  to  Step  5 

For  <j>  >  4,r 

Find  r  vs/2(l  -  sin  <£) 

x  =  r  cos  (X  -  X  ) 
o 

y  =  r  sin  (X  -  A.  ) 

(END) 

For  (^  <  <^T 

Initialize  E(x  -  XjT  ■  6y  =  0 
y  =  k  sin  3 

Find  A^(6X)  =  fiA(i  -  sin  (~<t>)) 

y  --  3  -  (X  -  X  ) 
o 

«2>  90°  +  y  -  3 

«y  =  {k  sin  fi  -  rT  •  sin  (0  -  y)}/N 
Use  SUBROUTINE  A 

to  find  (and  save)  0,  rT,  a,  b,  g,  c.  m,  p 
Find  Xj  -  X  +  6X 
Use  SUBROUTINE  A 

to  find  (and  save)  0^  rT1>  ay  gj,  c^  ni 
With  0,  ryi  a,  b,  g,  c,  m,  p,  y 
Use  SUBROUTINE  B 
to  find  (and  save)  x 

With  0j,  rT1,  a j,  bj.  gy  cy  my  py  y 
Use  SUBROUTINE  B 
to  find  (and  save)  x^ 
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Step  10. 


Find  (x  -  Xj)6v 


Kind  /  (x  -  x. ).  •  6v 


If  E  (x  -  x,  )<5v  <  A  (6 \)  go  to  Step  1 1 
1  •  0 

For  E(x  -  x.)6v  =»  A  ,(6v) 

1  0 


Head  (x,  y) 


Step  11.  For  £(x  -  x^)6y  <  A^(fiv) 
Set  y(new)  =  y(old)  -  6y 
Go  to  Step  8. 


(END) 


SUBROUTINE  A. 


Step  1. 


(Using  previous  entries,  including  A,  find  0,  r^,,  a,  b,  g,  c,  m,  p.  ) 
Initialize  a(u)  =  100 
a  (f)  =  0 


Find  7  =  3  -  X  +  A 

o 

0  =  90°  +  7-3 


*T(X)  -  (X  -  XQ) 


/ 2(1  -  sin  i 


»T(X)  -  0t(Y))/PO 


VT  =  k  sin  7 
VQ  =  k  cos  7  -  r^, 

A(,  =  27  (rad)  +  1/2  tan  yfVQ2  -  r 2  } 

7'  =  7  •  exp  (1/3  71/32  (3  -  7)} 
where  the  angles  in  the  brackets  are  in  radians. 
Step  2.  Increase  v  by  1. 

Find  a  =  (a(u)  +  a(l)}/2 

h  =/v/l  -  (1  -  VT/a) 


